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For the use of Gamma-Ray Bursts (GRBs) to probe cosmology in a cosmology- 
independent way, a new method has been proposed to obtain luminosity distances of 
GRBs by interpolating directly from the Hubble diagram of SNe la, and then calibrat- 
ing GRB relations at high redshift. In this paper, following the basic assumption in the 
interpolation method that objects at the same redshift should have the same luminosity 
distance, we propose another approach to calibrate GRB luminosity relations with cos- 
mographic fitting directly from SN la data. In cosmography, there is a well-known fitting 
formula which can reflect the Hubble relation between luminosity distance and redshift 
with cosmographic parameters which can be fitted from observation data. Using the Cos- 
mographic fitting results from the Union set of SNe la, we calibrate five GRB relations 
using GRB sample at z < 1.4 and deduce distance moduli of GRBs at 1.4 < z < 6.6 
by generalizing above calibrated relations at high redshift. Finally, we constrain the 
dark energy parameterization models of the Chevallier-Polarski-Linder (GPL) model, 
the Jassal-Bagla-Padmanabhan (JBP) model and the Alam model with GRB data at 
high redshift, as well as with the Cosmic Microwave Background radiation (CMB) and 
the baryonic acoustic oscillation (BAO) observations, and we find the ACDM model is 
consistent with the current data in 1-a confidence region. 

Keywords: Gamma rays: bursts; Cosmography 



1. Introduction 

Since an intrinsic relation between the peak luminosity and the shape of the light 
curve of SNe la has been found,^ SNe la has now been taken as near-ideal standard 
candles for measuring the geometry and dynamics of the universe. However, the 
maximum redshift of the SNe la which we can currently use is only about 1.7. On 
the other hand, the redshift of the last scattering surface of the cosmic microwave 
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background (CMB) is at z = 1091. 

Recently, Gamma-Ray Bursts (GRBs) were proposed to be a complementary 
probe to SNe la and CMB to explore the early universe. As the most intense ex- 
plosions observed in the uni vers e so far, GRBs are likely to occur in high-redshift 
range up to at least z — 8. 2. EE! Moreover, there are several luminosity relations of 
GRBs between the spectral and temporal properties which have been extensive dis- 
cussed, such as the isotropic energy (E'iso) - peak spectral energy (ii^poak) relatio n^ 



the luminosity (V) - spectral lag (riag) relation,^! the L - variability (y) relation,l^El 
the L - Epcak relation) ^'-*^*-*' the L - minimum rise time (trt) relation,^!!! and the 
coUimation-corrected energy {E^) - i?peak relation as well as several multiple re- 
lations such as the Eiso - Ep^ak - tb relation,^^ where ib is the break time of the 
optical afterglow light curves; the L - ii^pcak - 20.45 relation,l23] ^^ere T0.45 is the 
rest-frame "high-signal" timescale; and the L - Epcak - "Hag (or trt) relation.fl^ 
Many authors have made use of GRB luminosity indicators as standard 

candles at very high redshift beyond SNc la rcdshift range for cosmological 
research.IISIlMlDJ 12 20 21 23 24 13 26.27 2_8 29 30 31 32 33 34 35 36 37 38 39 

to the lack of the sample at low redshift which are cosmology independent, to 
calibrate the empirical GRB relations, one usually needs to assume a particular 
cosmological model with certain model parameters as a priori. As a result, the so- 
called circularity problem could prevent the direct use of GRBs for cosmology.l^Sl 
Many of works treat the circularity problem with statistical approach which car- 
ried out a simultaneo us fit of the parameters in the calibration curves and the 
cosmology! '*'^ ! ^^!^*-*!^'*' ^ '^^ ' However, the circularity problem can not be circumvented 
completely by means of the statistical approaches for an input cosmological model 
is still required in doing the joint fitting. 

More recently, Liang et al. proposed a new method to calibrate GRB luminos- 
ity relations in a cosmological mo del- independent way.^Sl The motivation of this 
calibration method is that objects at the same redshift should have the same lu- 
minosity distance in any cosmology. Thus the luminosity distance of a GRB at a 
given redshift can be obtained by interpolating directly from the Hubble diagram 
of SNe la, therefore GRB relations can be calibrated without assuming a particular 
cosmological model and the Hubble diagram of GRBs has been constructed. Follow- 
ing this cosmology-independent GRB calibration directly from SNe la, the derived 
GRB Hubble diagram can be us ed to constrain cosEaolo gical models at high red- 
shift avoiding circularity problem I44i 5 | 46 | 47 | 48 | 49 [50iaj. Capozziello & Izzo firstly 
used two GRB relations calibrated with the so-called Liang method to derive the re- 
lated cosmography parameters which related to the derivatives of the scale factor.ESl 
Liang et al. combined the updated distance moduli of GRBs obtained by the in- 
terpolating method with the joint data to find the contribution of GRBs to the 
joint cosmological constraints in the confidence regions of cosmological parameters, 
and reconstructed th e ac celeration history of the universe with the distance moduli 
of SNe la and GRBs.^^ On the other hand, besides the interpolation method, the 
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luminosity distance of a GRB can also be obtained directly from SNe la data by 
other mathematical approach. Liang & Zhang has proposed another approach to 
calibrate GRB relations by using an iterative procedure which is a non-parametric 
method in a model independent manner to reconstruct the luminosity distance at 
any redshift from SNe la.SS Similar to the interpolation method, Cardone et al. 
constructed an updated GRBs Hubble diagram calibrated by local regression from 
SNe la!^ Kodama et al. has proposed that the L — -Epoak relation can be cali- 
brated with one empirical formula fitted from the luminosity distance of SNe la!^ 
However, according to the formula fitting approach, various possible formula can be 
fitted from the SNe la data which could give different calibration results of GRB 
relations. As the cosmological constraints from GRBs are sensitive to GRBs cal- 
ibration results, and the fitting procedure depends seriously on the choice of the 
formula, the reliability of this method should be tested carefully. In other words, 
we should find one certain formula which is totally independent of any cosmological 
models and could accurately evaluate the Hubble relation. 

In Cosmography,"^ there is a well-known formula reflecting the Hubble rela- 
tion between luminosity distance and redshift which can be extracted directly from 
basic cosmological principles and observation data, with cosmography parameters 
(the deceleration, jerk and snap parameters: q, j, and s) which are only related to 
the derivatives of the scale factor without any priori assumption on the underlying 
cosmological model. Recently, several authors have already used the cosmographic 
parameters fitting from SNe la and/or GRBs dataset to constrain cosmological 
parameters .f^"' If viewing this point from another angle, the cosmographic 

formula can be considered as a perfect fitting function to calibrate the GRB rela- 
tions using SNe la data, as long as we take the same assumption that objects at 
the same redshift should have the same luminosity distance in any cosmology. In 
this paper, instead of the interpolation method using in Ref. 1431 we propose another 
new approach to calibrate GRB luminosity relations with cosmographic fitting from 
SNe la data. The structure of this paper is arranged as follows. In section 2 we give 
a brief review of the cosmographic Hubble relation between luminosity distance and 
redshift. In section 3, we calibrate five GRB luminosity relations with cosmographic 
fitting results from SNe la data. In section 4, we construct the Hubble diagram 
of GRBs obtained by using the cosmographic methods and constrain the da rk en- 
ergy parameterization models of the Chevallier-Polarski-Linder (CPLjmodel,^^^^^ 
the Jassal-Bagla-Padmanabhan ( JBP) model^^ and the Alam model^ with GRB 
data at high redshift, as well as with the Cosmic Microwave Background radiation 
(CMB) and the baryonic acoustic oscillation (BAG) observations. Conclusions and 
discussions are given in section 5. 
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2. Cosmographic Hubble relation between luminosity distance and 
redshift 

A completely new cosmology branch ~ cosmography, in which framework cosmology 
is pure kinematics and completely independent of the underlying dynamics gov- 
erning the evolution of the universe has been introduced since 1970s!^ The only 
assumption is the basic symmetry principles (the cosmological principle) that the 
universe can be described by the Friedmann-Robertson- Walker metric. By means 
of a Taylor series expansion of the luminosity distance, Visser gave a formulation 
that the luminosity distance can be expressed as a power series in the redshift up 
to the forth order;"^^ 
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where the coefficients of the expansion are the so-called cosmographic parameters, 
Hubble parameters {H = a{t)/a{t)), deceleration parameters {q = —H~^a{t)/a{t)), 
jerk parameters {j = H^'^a^'^\t)/a{t)), and snap parameters (s = H^'^a^^^ {t)/a{t)), 
which related to the scale factor a{t) and its higher order derivatives (the subscript 
"0" indicates the present value of the parameters) . Obviously pure cosmography by 
itself will not predict anything about the scale factor a{t), we have to turn to the 
observational data such as SNe la to infer the history of the scale factor a{t) and 
some important information about expanding history of our universe. 

In order to avoid problems with the convergence of the series for the highest 
redshift objects as well as to control properly the approximation induced by trun- 
cations of the expansions, Cattoen & Visser pointed out that it is useful to recast 
as a function of an improved parameter y = z/{l + z) and constrained the cosmo- 
graphic parameters using SNe la data.l^In such a way, being z e (0, oo) mapped 
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where flo = 1 + / H^a?{to) is the total energy density. For the flat universe, Vlo = 
1. The luminosity distance as the logarithmic Hubble relations can be expressed as: 
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therefore the distance modulus can be given by 
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Recently, Vitagliano et al. fitted two different truncations (Cosmography I: with- 
out the third order term {y^); Cosmography II: with the third order term {y^)) of 
the above expansion with the SNe la and GRB datasetsfS^ With a flat universe 
(rio = 1) prior, the fitting results showed that the Union dataset!^ gave more strin- 
gent constraints on the parameters. For Cosmography I, the cosmographic fitting 
results with the Union dataset are 

go = -0.58 ±0.24 , jo + f^o = 0.91 ± 2.21 , 

and for Cosmography II, the cosmographic fitting results with the Union dataset 
are 

go = -0.50 ±0.55 , jo + ^^0 = -0.26 ±9.0, sq = -4.13 ± 129.79 . 

It is noted that the Union dataset of 307 SNe la didn't include the 90 SNe 
la data from CfA^^ due to their extremely low redshift (z < 0.1), which would 
not affect the calibrated results for GRB luminosity relations at z > 0.17.^^ With 
the above cosmographic fitting results from the Union dataset, the cosmographic 
Hubble relation can be considered as a perfect function to deduce the distance 
moduli of GRBs directly from SNe la data. In the next section, we will deduce the 
distance moduli of GRBs and then calibrate the GRB luminosity relations with 
Cosmography I and II respectively. 

3. The Calibration of the Luminosity Relations of Gamma-Ray 
Bursts 

We adopt the 69 GRBs provided in Ref. [30l as our sample for calibrating the GRB 
luminosity /energy relations. We first deduce the distance moduli of GRBs at z < 1.4 
within our sample. Then using these deduced distance moduli and the redshifts of 
corresponding GRBs, we calibrate five GRB luminosity/energy relations i.e., the riag 
- L relation, the V - L relation, the L - £'peak relation, the - -Epeak relation and 
the trt - L relation. These luminosity relations of GRBs can be generally written 
in the form 

logy = a + 61oga;, (5) 

where a and h are the intercept and slope of the relation respectively, in addition we 
introduce c as the linear correlation coefficient of the relation which will be calcu- 
lated together with a and b below; y is the luminosity (L/ergs~^ or energy E-y/eig); 
X is the GRB parameters measured in the rest frame, e.g., Tiag(l + z)^^/(0.1 s), 
V{1 + z)/0.02, i;pcak(l + z)/(300 keV), rRT(l + 2)^7(0.1 s), for the corresponding 
relations above. For the x values, we adopt the data from Ref. '30| for the y val- 
ues, we drive them with the adjusted luminosity distance of GRBs calculated with 
cosmographic fitting method (Cosmography I and II). The isotropic luminosity of 
a burst is calculated by 

L = 47r4Pboio, (6) 



4, 2012 1:26 WSPC/INSTRUCTION FILE Cosmography 



6 

where Pboio is the bolometric flux of gamma-rays in the burst and cJl is the lu- 
minosity distance of the burst. The isotropic energy released from a burst is given 

by 

Si.o = 47r45boio(l + (7) 

where S'boio is the bolometric fluence of gamma-rays in the burst at redshift z. The 
total coUimation-corrected energy is then calculated by 

where the beaming factor, Fbcam = (1 ~ cos^jct); and the value of the jet opening 
angle fljet is related to the jet break time (ib) and the isotropic energy for an 
Earth-facing jet, i?7,iso,52 = £'7,iso/10'^^erg. When calculating E'^^so, we also use 
the cosmographic fitting method from SNe la to avoid the circularity problem. 

We determined the values of the intercept (a) and the slope (6) with their \-u 
uncertainties with the same regression method (the bisector of the two ordinary 
least-squares) used in Ref. |30l |43l The calibrate results for Cosmography I and 
II are summarized in Table 1. From Table 1, we find that the calibration results 
obtained using two cosmographic fitting (Cosmography I and II) are fully consistent 
with each other. We also find that results obtained by the cosmographic methods 
differ only slightly from, but still fully consistent with those calibrated by using the 
interpolation method and the iterative procedure with the same GRB sample (see 
details inthere^^ ^4). 

Table 1. Calibration results (for a=intercept, 6=slope, c=correlation coefficient) with 
their l-c uncertainties, for the five GRB luminosity /energy relations within the sample 
at z < 1.4, using two cosmographic fitting results (Cosmography I and II) directly from 
SNe la data. 





Cosmography I 




1 Cosmography II 




Relation 


a 


b 


c 


a 


b 


c 




52.13±0.10 


-I.10i0.13 


-0.89 


52.21±0.11 


-1.13±0.14 


-0.88 


V-L 


52.47±0.13 


2.02±0.27 


0.64 


52.54±0.13 


2.06±0.27 


0.65 




52.14±0.09 


1.64±0.10 


0.89 


52.21±0.09 


1.67±0.10 


0.89 




50.83±0.06 


1.87±0.11 


0.95 


50.88±0.07 


1.93±0.10 


0.95 


trt-L 


52.51±0.11 


-1.29±0.11 


-0.77 


52.58±0.11 


-1.31±0.12 


-0.77 



4. The Hubble diagram of gamma-ray bursts 

With the cosmographic fitting results (Cosmography I or II), the moduli for the 27 
GRBs at z < 1.4 can be directly obtained from SNe data, therefore we calibrate 
GRB luminosity relations in a completely cosmology-independent way. Further- 
more, if assuming that GRB luminosity relations do not evolve with redshift, we 
are able to obtain the luminosity (L) or energy (E^) of each burst at high redshift 
(z > 1.4) by utilizing the calibrated results from Cosmography. Consequently, the 
corresponding luminosity distance (cIl) can be derived from Eq.(12) ^ Eq.(14) and 
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the corresponding distance modulus can be calculated as /i = 5 log(dL /Mpc) + 25. 
The uncertainty of the value of the luminosity or energy deduced from a GRB 
relation is 

^logy - + (^6 ^ogxf + (OASma^/xf + a^,, (9) 

where CTo, (T& and ax are l-cr uncertainty of the intercept, the slope and the GRB 
measurable parameters, and agys is the systematic error in the fitting that accounts 
for the extra scatter of the luminosity relations!^ Note that the uncertainty of 
modulus for each luminosity indicator depends on whether Pboio or S'boio is used: 

= [(2.5aiogL)' + (1.086ap,„,jPboio)']'/', (10) 

or 

= [(2.5(7iogi;j' + (1.086as,„,J5bolo)' + (1.086(7;^,_7i^^beam) Y^'- (11) 

For five luminosity indicators, each burst will have up to five estimated distance 
moduli, we hence use the same method in Ref. [301 to obtain the best estimated fi 
for each GRB which is the weighted average of all available distance moduli: 

A'=(EA'i/0/(E^;.')' (12) 

i i 

with its uncertainty = (^^ "'^7;^)^^^^: where the summations run from 1 to 5 over 
the five relations used in Ref. 1301 with available data. Until now we have ultimately 
obtained the 42 GRB moduli at z > 1.4 by utilizing the five relations calibrated 




0123456789 



Fig. 1. Hubble Diagram of 307 SNe la (red dots) and the 69 GRBs (circles) obtained using the 
cosmographic method. The 27 GRBs at 2 < 1.4 are obtained by the cosmographic fitting from 
SNe la data {black circles), and the 42 GRBs at z > 1.4 {blue circles) are obtained with the 
five relations calibrated with the sample at z < 1.4 using the cosmographic method {blue circles: 
Cosmographic I; blue stars: Cosmographic II). The curve is the theoretical distance modulus in 
the concordance model {w = —1, Omo = 0.27, Ha = 0.73), and the vertical dotted line represents 
z = 1.4. 
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with the sample at z < 1.4 using the fitting method. We have plotted the Hubble 
diagram of 307 SNe la and the 69 GRBs obtained using the cosmographic methods 
in Fig 1. 

Constraints from the GRB data can be obtained consequently by fitting the dis- 
tance moduli fJ.{z). The value of the observed distance moduli can be calculated 

by 

^t^^l^ -2 ' (13) 

where /iobs(^i) is the observed distance modulus for the GRBs at redshift Zi with 
its error cr^j; n{zi) is the theoretical value of distance modulus from a dark energy 
model. As mentioned above, high-redshift GRBs are rare database for constraining 
cosmological parameters. Since those 42 GRBs' moduli (1.4 < z < 6.6) calculated 
above are completely cosmological model independent, we thus utilize this dataset 
together with CMB and BAO observations to constrain specified cosmological mod- 
els. 

For the CMB observation, we choose one shift parameter R to limit the model 
parameters. In a fiat universe, it can be expressed as 



where the last scattering redshift zis — 1091.3 from the 7-year WMAP results, and 
the observational value R — 1.725 ± 0.01&^. The Xcmb value is 

_ (i^- 1.725)^ 

XCMB — Q g^g2 • 

For the BAO observation, the size of bary on acous tic oscillation peak can be 
used to constrain the cosmological parameters.l^SIMEZl This peak can be denoted 
by a parameter A, which can be expressed aJ^ 



^BAO Jo 



E{z 



where zbao = 0.35. The observational value is A = 0.469(ns/0.98)"" -^5 ± 0.017,1^ 
with the scalar spectral index rig = 0.963 from WMAP7 data.l^l The Xbao value is 

2 _ [A~ 0.469(n,/0.98)-°-35]2 
Xbao- OOIT^ ' ^ ' 

Here we combine these two probes with the GRBs dataset above by multiplying 
the likelihood functions. The total value is 



Xtotal — XgRB + XCMB + XbAO- 

(18) 
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In the following, we will show the constraining progress and results for three 
dark energy parametrization models: CPL model, JBP model and Alam model. 
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Fig. 2. Contours of likelihood in the (uio,iOa) plane in the CPL dark energy model for a flat 
universe, from CMB and BAO observations together with the 42 GRBs data (z > 1.4) obtained 
by the cosmographic method. The red line is for Cosmography I, and the red plus sign denotes 
the best-fit values (uiq = —0.79, Wa = —0.36). The black line is for Cosmography II, and the black 
plus sign denotes best-fit values (wq = —0.91, Wa = —0.06). The contours correspond to 1 and 2-a 
confidence regions. 

(1) Chevallier-Polarski-Linder (CPL) model,^^^^^ in which dark energy with a 
parametrization EoS: 

w{z)^Wo+Wa— . (19) 

1 + Z 

The corresponding luminosity distance for a flat universe is 

fly 

dL=cH^\l + z) (20) 



Table 2. The best-fit values of wq and Wa (or A\ and A2) with l-cr uncertainties, as well as the best-fit values of f2ji,/, 
Hq, and Xmin' CPL model, JBP model and Alam model, with GRB data (calibrated with Cosmography I and II 

fitting from SN la data), as well as with CMB and BAO observations. 





CPL Model 


,IBP Model 


Alam 


Model 




Cosmography I 


Cosmography II 


Cosmography I 


Cosmography II 


Cosmography I 


Cosmography II 


wo{Ai) 


' "-0.34 


-0 9l+"-^« 
-0 06+^ "" 


-0.93tl\ll 


-i.37li!:- 

^■"^^-3.48 


60+" "'^ 
"■""-0.97 

-0.llto.25 


-0.39l-.«'; 

o.iotS:?! 




0.29 


0.28 


0.29 


0.25 


0.29 


0.26 


Ho 


80 


75 


80 


79 


79 


78 




28.03 


26.91 


28.02 


26.91 


28.03 


26.90 




Cosmography I 
+ Best fit for Cosmograptiy I 
■',':V;H'--Cosmograptiy II 
+ Best fit for Cosmograptiy 1 1 
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where 

E{z) - [(1 + zf^M + (1 - f^Af)(l + z)3(l+-o+-a)g-3«,„./(l+^)] (21) 

We set Hq and 51^/ as free parameters and their prior values are 
64km s"^Mpc"^ < Hq < 80km s^^Mpc"^ and 0.23 < I^a/ < 0.31 respectively. 
We also assume this prior in the following analysis. For cosmography I, we find 
Xmin = 28.03 at {Ho , rJA/}={80km s"^Mpc"\ 0.29}, the best fitting values 
are Wq = — 0.791q34 and Wa = —0.36~^\y^ (If'')- For cosmography II, we find 
xl^i^ = 26.91 at {Ho , flM}={75km s"^Mpc"\ 0.28}, the best fitting values are 
Wo = — 0.91^q'32 and Wa = —0.06t_\'QQ (If). Fig. 2 shows the constraints on wo and 
Wa parameters for the CPL model. We present the best-fit value of wo and Wa with 
la uncertainties, as well as the best- fit values of ^m, Ho, and Xmiw ^o^' ^he CPL 
model in Table 2. 

(2) Recently, Jassal, Bagla & Padmanabhan argued that in CPL model some 
problems will present at high redshifts, they thus proposed a new paramctrization 
EoS (JBP model) ^ 

Z 

w{z) ^ Wo + Waj—^ (22) 

(1 + zy 

which can model a dark energy component that has the same value at lower and 
higher redshifts, with rapid variation at low z. The corresponding luminosity dis- 
tance for a flat universe is the same with Eq.(26), and Eq.(27) becomes 



E{z)= [(l + z)-''r!Af + (l-f^Af)(l + z)'^'+'""^e3'"»"'/2(i+.)Y/2^ (23) 

For cosmography I, we find Xmin = 28.02 at {Ho , 17A/}={80km s"^Mpc-\ 
0.29}, the best fitting values are wo = -0.74^°;^^ and Wa = -0.93t^;5^ (la). For 
cosmography II, we find Xmin = 26.91 at {Ho , f7Af}={79km s-^Mpc-\ 0.25}, the 
best fitting values are wo = — 1.37to.5i and Wa = 2.02t3;4g (la). Fig. 3 shows the 
constraints on wo and Wq parameters for the JBP model. We present the best-fit 
value of Wo and Wa with la uncertainties, as well as the best-fit values of 51 a/, Ho-, 
and Y^:„, for the JBP model in Table 2. 

(3) The third dark energy paramctrization model that we consider is^ 

MA = l±L^l±2^^ _ 1, ,24) 

where 

^Ide(^) = ^i(1 + z) + ^2(1 + z)2 + 1 - - - ^2, (25) 

thus 

E{z) = [(1 + zf^M + Ai(l + z) + ^2(1 + z)2 + 1 - - - (26) 

For cosmography I, we find Xmin — 28.03 at {Ho , 57A/}={79km s^^'^Mpc^^, 
0.29}, the best fitting values are Ai = 0.6lo;g7 and yl2 = -O.lllo.l (1^^). For 
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cosmography II, we find xf„i„ = 26.90 at {Hq , fJM}={78km s"^Mpc"\ 0.26}, the 
best fitting values are Ai = -OM^oIt and A2 = O.llo;2l(lcr). Fig. 4 shows the 
constraints on Ai and A2 parameters for the Alam model. We present the best-fit 
value of Ai and A2 with la uncertainties, as well as the best-fit values of 51^/, Hq, 
and Xmin' for the Alam model in Table 2. 

We can find that when combining the GRBs dataset obtained by cosmographic 
fitting method with the CMB and BAO observations, the combined data can give 
more stringent results on different dark energy parametrization models. It is noted 
that all these results for above models are consistent with the ACDM model in 1-a 
confidence region. 

5. Summary and Discussion 

Due to the lack of the GRB sample at low redshift, there has been a so-called 
circularity problem which can always be a obstacle for applying GRBs data to 
constrain cosmological parameters. Based on the basic assumption that objects at 
the same redshift should have the same luminosity distance, Liang et al. proposed 
a new method to calibrate GRB relations in a completely cosmology-independent 
way, namely obtaining the distance modulus of a GRB by interpolating from the 
Hubble diagram of SNe la and then calibrate the GRB relations with these cal- 
culated distance modulPSl. There is a well-known fitting formula in cosmography, 
which can reflect the Hubble relation between luminosity distance and redshift with 
cosmographic parameters which can be fitted from SNe la. In this work, we propose 
another approach to calibrate GRB luminosity relations with cosmography fitting 
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Fig. 3. Same as Fig. 2, but fitting the JBP model. For Cosmography I: the best-fit values {wq = 
— 0.74, = —0.93); for Cosmography II: the best-fit values (uiq = — 1-37, = 2.02). 
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from SNe la data. We adopted the fitting results from the Union set of SNe la for 
the so-called Cosmography I and Ilf^and calibrate five GRB relations by this cos- 
mographic fitting method. The calibration results obtained using two cosmographic 
fitting are fully consistent with each other. Assuming that GRB luminosity relations 
do not evolve with redshift, we obtained the distance modulus of the GRB data at 
higher redshift 1.4 < z < 6.6. Note that the circularity problem could be completely 
avoided if we apply these GRBs data to constrain cosmological parameters. We thus 
constrained three parameterization dark energy models of the GPL, JBP and Alam 
models by combining the GRB data at high redshift with the GMB and BAO ob- 
servations, and the constraint results are all consistent with the AGDM model in 
1-cr confidence region. 

Compared to previous cosmology-independent calibration method, we find that 
results obtained by the cosmographic method consistent with those calibrated by 
using the interpol ation m ethod and the iterative procedure from SNe la with the 
same GRB sampl(S2lM 

this situation are consistent and similar with the recent 
conclusions that the cosmographic Amati relation agrees in the errors with other 
cosmology-independent calibrations,'^ and that the calibration coefficients and the 
intrinsic scatter actually do not depend on the adopted calibration procedure se- 
riously, based on the use of a fiducial model, the cosmographic method and the 
local regression.!^ It should be noted that the fitting procedure used in Kodama 
et al. 1^ depends seriously on the choice of the formula. In Cosmography, the cos- 
mographic formula is totally independent of any cosmological models and could 




Fig. 4. Same as Fig. 2, but in the {Ai,A2) plane in the Alam model for a flat universe. For 
Cosmography I: the best-fit values {Ai = 0.6, A2 = —0.11); for Cosmography II: the best-fit 
values (Al = -0.39, A2 = 0.1). 
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accurately evaluate the Hubble relation between luminosity distance and redshift, 
and can be considered as a perfect fitting function to calibrate the GRB relations 
using SNe la data. The reliability of this method should be more reasonable than 
the fitting procedure which chooses a formula in arbitrary!^ However, since the dis- 
tance moduli of GRBs calculated by the fitting method including the cosmographic 
method are obviously relevant to all of SN la data points for calibration, it will be 
not appropriate fo r direc tly combining these two datasets to constrain the cosmo- 
logical parametersP^^^ But like showing in this work, we can combine the GRBs 
dataset obtained by cosmographic fitting method with the CMB and BAO observa- 
tions an d th e combined data can give more stringent results. As to the interpolation 
method,S2l for deducing one individual GRB, only a few SNe la close to this GRB 
were used. In this case, the measurement error for a single SNe la's modulus may 
infiuence the final result of calibration. In the cosmographic method, the full infor- 
mation of SNe la data is completely used in the calibration. It should be noted that 
our approach provides no more accurate cosmological parameter constraints than 
other works such as the simultaneous fit method. However, the primary motiva- 
tion of this work is not on improvement of the statistical error, rather on avoiding 
the circularity problem more clearly in logic. 

For comparing GRB to SNe la, GRBs are almost immune to dust extinction, 
whereas SNe la observations suffer extinction from the interstellar medium. GRBs 
can extend the Hubble diagram to much higher redshifts beyond SNe la data. Differ- 
ent dark energy models may have very different Hubble diagrams at high redshifts. 
On the other hand, due to the large statistical scatters of the relations and the 
small dataset compared to SNe la, the contribution of GRBs to the cosmological 
constraints would not be sufficiently significant at present. However, it should be 
noted that a single GRB at high redshift will provide more information than a sin- 
gle maximal redshift SNe la,"^ and the larger scatter of GRB data is somewhat 
coming from the observational limit of nowadays technological level, which means 
the GRB data may eventually turn into good data as our observational capability 
enhance. It has been found that cosmological constraints would improve substan- 
tially with more simulated GRBs expected by future observations through Monte 
Carlo simulations.!^ With more and more GRBs observed from Fermi Gamma-ray 
Space Telescope with much smaller scatters, and its combination with the increasing 
Swift data, GRBs could be used as an additional choice to set tighter constraints 
on cosmological parameters of dark energy models. At that time when GRB data 
becomes common used data sample to constrain high redshift cosmology, the so 
called circularity problem for GRB data will be a key problem to solve. Moreover, 
it is noted that accurately calibrating the GRB relations is also very important for 
the GRB theory system itself even without considering the cosmology. We have to 
note stress again that the high-redshift GRB dataset we obtained here is completely 
cosmology-independent, it will ultimately fill the data crack between SNe la and 
CMB after more GRBs being observed in the future. 
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